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Abstract 

We present new deep K' counts of stars at the South Galactic Pole (SGP) taken 
with the NAOJ PICNIC camera to if' =17.25. Star-galaxy separation to iC=17.5 
was accomplished effectively using image profiles because the pixel size we used is 
0.509 arcsec. We interpret these counts using the SKY (Cohen 1994) model of the 
Galactic point source sky and determine the relative normalization of halo-to-disk 
populations, and the location of the Sun relative to the Galactic plane, within the 
context of this model. The observed star counts constrain these parameters to be: 
halo/disk ~l/900 and z Q =16.5±2.5 pc. These values have been used to correct 
our SGP galaxy counts for contamination by the point source Galactic foreground. 
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1. Introduction 

The advent of new near-infrared (NIR) cameras is providing a wealth of valuable star 
counts that can be used to constrain models of the Galaxy. These models are critical 
to the quest for any near- and mid-infrared cosmic background radiation. Recently 
Minezaki et al. (1997: hereafter designated as Paper I) have presented differential K'- 
band galaxy counts very close to the South Galactic Pole (SGP) taken with the PICNIC 
camera using a pixel size of 0.509 arcsec on the the Australian National University's 2.3 m 
telescope at Siding Spring Observatories, Australia. The purpose of the present paper is 
twofold: to describe our efforts to interpret and model the stellar content of these counts 
through use of the SKY model (Cohen 1994; Wainscoat et al. 1992, hereafter referred 
to as WCVWS) and to justify the manner in which we removed the stellar foreground 
contribution to these galaxy counts in Paper I, to isolate the extragalactic component. 

We very briefly summarize the relevant observational details of Paper I's "bright sur- 
vey" in Section 2, and the version of the SKY model that we use in Section 3. In Section 
4 we discuss the way in which we have separated the interpretation of the star counts 
into the measurement of the halo/disk normalization factor and of the location of the 
Sun with respect to the Galactic plane. In this we follow precisely the approach taken by 
Cohen (1995), who pursued the same two parameters in the context of SKY, primarily 
through the use of IRAS source counts at the two Galactic Poles and of high latitude 
Shuttle-based far-ultraviolet star counts. Our values for these parameters are in good 
agreement both with Cohen's values and those found by conventional techniques and 
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modeling in the visible. Section 5 presents an analysis of our estimated uncertainties in 
the extrapolation of star counts to faint K' magnitudes using SKY. 

2. The observations 

A full description of the observations, and other analyses, were presented in Paper I. The 
"bright survey" was carried out during August and September 1994, using the Australian 
National University's 2.3 m telescope at Siding Spring Observatory, Australia, equipped 
with the PICNIC near-infrared camera (Kobayashi et al. 1994), developed at National 
Astronomical Observatory, Japan. PICNIC uses a NICMOS3 array (256 x 256 pixels) 
with a field of view of 2.2 x 2.2 arcmin 2 and a pixel scale of 0.509 arcsec pixel -1 . In 
order to reduce the thermal sky background, we used a K' filter, which has the same 
transmission curve as the 2MASS K s filter (McLeod et al. 1995). 

The survey was centered at B2000 coordinates {a, 5) = (0 h 50 m 48 s , -27°43'34") or 
(l,b) = (316.27°, -89.39°), and covered 180.8 arcmin 2 . The FWHM of PSF was 1.5 
arcsec and the limiting magnitude with 80% completeness was determined as K' — 19.1 
from the simulations of the recovery of stellar profiles (cf. Paper I). 

Stars and galaxies in the survey were separated to 1.5 mag brighter than the limiting 
magnitude, K' = 17.5, based on two morphological parameters, the FWHM and the 
irl, where the FWHM was measured by Gaussian fitting of the radial profile by IRAF 
imexamine task and the irl was the intensity-weighted first moment radius which was 
measured by FOCAS. 

The K' band star counts to K' = 17.25 obtained from the survey are tabulated in 
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Table 1. 



3. The version of SKY used 

We have utilized version 4 of the SKY model, exactly as documented by Cohen (1994) 
and Cohen et al. (1994). The difference between implementations of the "standard" 
K and the newer K' filters in SKY have been explored by Cohen (1997), who found 
differences of at most 0.03 mag between the absolute magnitudes, and M K i, for 
the entire set of 87 categories of source with which SKY populates its representation of 
the sky. These differences are negligible compared with both the intrinsic dispersions 
and uncertainties in the absolute magnitudes of the many types of celestial source, and 
with the photometric uncertainties in the measurement of K' at faint magnitude levels. 
Therefore, our simulations with SKY were based on SKY's "hardwired" standard K 
passband. 

SKY incorporates five geometric components: the (thin) disk; bulge; spiral arms, 
local spurs, and Gould's Belt; molecular ring; and halo. SKY does not formally include 
a thick disk (cf. Yoshii 1992; Yoshii et al. 1988), principally because it has negligible 
influence on currently available infrared star counts all of whose Poisson errors greatly 
exceed the thick disk's expected contribution in any direction. 

4. The SKY simulations 

We first motivate the analysis with SKY by a simple analytical approach. Assuming 
that the vertical disk scale-height h z is the same above and below the disk plane, the 
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solar displacement z can be roughly evaluated from the ratio between the observed star 

counts in the north and south Galactic poles (hereafter referred to as NGP and SGP): 

NGP e- x °(x 2 Q + 2x Q + 2) 
SGP ~ 4-e- x °(xl + 2x Q + 2) ' 

where x & = z Q /h z . We currently lack comparably deep K' counts at the NGP to combine 

with our own SGP infrared star counts. The only counts we have examined near the 

NGP in K' are from the 2MASS pilot study of SA57 at 6=85° (Skrutskie 1996), to a 

limit of 14th magnitude. Using our own and these 2MASS data we can determine only 

that z is greater than zero, because so few stars are available to the common limit of 

14. 

Therefore, we attempted to derive x Q and, hence, z Q (since h z is known for each 
category of source) using the relation 

POO 

SGP = up h 3 z / e~ x x 2 dx = up Q hl{A - e~ x ° {x% + 2x Q + 2)) , 

Jxq 

for one-sided star counts, to constrain z Q since p is known, uj represents the area in 
steradians over which we assess the star counts, and we summed over all categories of 
star. 

Using values for scale height and space density purely for the thin disk from WCVWS 
(their Table 2), and our observed cumulative counts (corrected for incompleteness) of 
67.17 stars within the bright survey area, we were again unable to obtain more than an 
approximate value for z Q because of the appreciable Poisson uncertainties in our counts. 

Consequently, we chose to interpret the counts strictly within the context of SKY in 
which many stellar populations are represented, along with the pertinent geometry of 
halo and disk components. We compared the counts predicted by SKY with the star 
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counts in Table 1, after applying the corrections for incompleteness (hereafter referred to 
as "corrected counts"), and after removal of galaxies on the basis of their much broader 
photometric profiles. These counts were binned into 0.5 mag wide bins, then converted 
to differential counts per square degree per magnitude interval. The sensitivity of our 
analysis in determining the optimal value of z Q is less than that associated with halo/disk 
so we first sought the best solution for halo/disk. We fixed the Sun's position at 15 pc 
north of the Galactic plane (cf. Cohen 1995) and ran a series of predictions for differential 
counts specific to the i-T-band at the central coordinates of the mosaic of PICNIC fields, 
varying halo/disk from to 2, in steps of 0.1, where the unit is the standard ratio of 
1/500 adopted by WCVWS. Following Cohen (1995), we analyzed three statistics to 
select the best solution for halo/disk and z Q , namely: unweighted algebraic and inverse- 
variance-weighted average deviations (AD and WAD, respectively) and the unweighted 
sum-squared-deviations (SSD). The former require the derivation of zero crossings; the 
latter of an absolute minimum. The three distinct quantities were interpolated onto a 
much finer grid of halo/disk ratio. The formal x 2 analysis was again (cf. Cohen 1995) 
deemed to be unsuitable here because the dynamic range of the observations is close to 
a factor of 100 in both axes, leading to inappropriately high weighting of the faintest 
points observed. Our best fitting values were: WAD, 0.56; AD, 0.59; SSD, 0.54 (in units 
of the standard ratio of 1/500 used by WCVWS). 

These procedures lack a natural associated error. To extract the formal uncertainties 
we assessed the sensitivity of this technique by increasing the corrected counts by 1 star 
in the faintest 3 (0.5 mag wide) bins, and similarly decreasing those in the brightest 3 
bins by 1 star. We then repeated the procedure, switching the senses of the increase 

7 



and decrease. (Randomly adding or subtracting stars would lead only to a best-fitting 
set of counts with more or less the same shape as that observed, hence the same value 
of the parameter would result.) In this way we derived plausible estimates of the non- 
systematic errors associated with the three statistics as: WAD, ±0.015; AD, ±0.018; 
SSD, ±0.031 (strictly, these errors are slightly asymmetric in the positive and negative 
directions). We, therefore, derive the optimal value for halo/disk as 0.56±0.03 in units 
of the standard ratio of 1/500. 

We fixed the ratio halo/disk at this formal best-fitting value and similarly constructed 
the three statistics for z by running a series of SKY predictions differing by 1 pc in their 
input value of z , covering the range 5-30 pc. The identical methodology, applied to the 
analysis of z , yielded: WAD, 14.6; AD, 18.9; SSD, 16.0 pc, with uncertainties (assessed 
in the manner described above) of ±3.6, 2.0, and 1.7 pc, respectively. Therefore, we take 
16.5±2.5 pc as the best estimate for the solar location using the PICNIC counts and 
SKY version 4. 

Figure 1 compares this best-fitting simulation with our corrected star counts. This 
simulation was extrapolated to assess the contribution of faint stars to the total counts 
in the PICNIC "bright survey" at K' >17.5, below which we did not attempt to separate 
stars and galaxies. 

5. Uncertainties in extrapolating counts with SKY 

We have attempted to quantify the likely uncertainties in SKY's star counts as follows. 
There are three primary contributors to these uncertainties; two are traceable directly 
to the separate influences of the imprecision with which we know z and halo/disk. The 
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third relates to our knowledge of the space density of the dominant contributors to the 
total differential counts at the faintest K' magnitudes. 

To assess the first two elements of uncertainty, we ran SKY for the four extreme 
cases in which one of z and halo/disk was held at its optimal value (Section 4), while 
the other was set to its formal mean±3<7 level. Consequently, these four predictions 
were for (z , halo/disk) of (16.5,0.47), (16.5,0.65), (9,0.56), and (24,0.56). From the half- 
ranges of the difference in total counts we derived the separate fractional uncertainty 
contributions from the two parameters. Figure 2 includes these contributions over the 
relevant K' range. 

SKY returns a highly detailed "log" file of its predictions that breaks down the total 
differential counts into the five geometric components, and each of those into the distinct 
counts from all the 87 categories of source over any desired range of magnitudes. From 
these files, and on the basis of simulations run against other sets of observed counts 
in high latitude fields kindly supplied by Hammersley (1996: the North Galactic Pole 
complete to H ~15) and by Skrutskie (1996: for 6=85° from 2MASS in JHK\ complete 
to ~16, 15, and 14, respectively), as well as those obtained by Meadows (1994: complete 
to K' ~17.5), (i) we have determined that the "M LATE V" stars fulfil the role of 
this dominant population, increasingly so as one goes to fainter K' magnitudes; and (ii) 
we estimate the likely uncertainty in space density of these stars at ±0.2 dex. We can, 
therefore, derive the impact upon total counts of a ±0.2 dex change in the space densities 
of these cool dwarf stars alone. This provides the third component of uncertainty in our 
analysis (see Figure 2). 

All three components are distinct and independent so we combined them, as fractional 
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uncertainties in a root-sum-square fashion, to yield our final estimate (the curve labeled 
"Total" in Figure 2) of the ±3<r fractional uncertainty associated with the extrapolation 
of the PICNIC corrected counts to faint magnitudes using SKY's predictions. 

By extrapolating the fitted SKY model to K' >17.5, the contributions of the star 
counts to the total counts were estimated as about 7% at K'=18, and 5% at K'=19, 
which were just comparable to the Poisson errors of galaxy counts in the PICNIC "bright 
survey". Since the estimated uncertainty of SKY predictions is about 0.16 dex, or 40%, 
at K' =18-19, the uncertainties of the star count predictions would not affect the galaxy 
counts of the PICNIC "bright survey" presented in Paper I. 

6. Discussion 

There is substantial corroborative evidence in favor of an offset of the Sun from the 
plane of order 20 pc from quite different disciplines (see the references in Cohen 1995), 
such as the distributions of optically-known Wolf-Rayet stars (Conti & Vacca 1990: 15 
pc) and of diffuse Galactic infrared radiation (Hammersley et al. 1994: 15.5 pc; Arendt 
et al. 1996: 18 pc), all of which are somewhat smaller than the value found from analysis 
of faint red stars (Yamagata & Yoshii 1992: 40 pc). We conclude that our best value, 
from these deep near-infrared counts within the context of the SKY model, is about 16.5 
pc, with a formal la uncertainty of about 2.5 pc. 

WCVWS originally adopted a halo/disk normalization factor of 1/500 from Bahcall 
& Soneira (1984), within the estimated uncertainties obtained by Schmidt (1975), who 
derived a best value of 1/800 with rough probable range 1/550 to 1/1200. Bahcall et al. 
(1983) similarly deduced ratios between 1/200 and 1/1500, and Cohen (1995) suggested 
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a number density ratio of halo/disk of 1/1250 for SKY version 4, at the lower bound of 
Schmidt's suggested approximate range, but within that determined by Bahcall et al. 
(1983). The uncertainty inherent in Cohen's (1995) estimate of halo/disk is probably 
at least ±0.15, in units of WCVWS's value (1/500). Within the context of SKY, the 
PICNIC counts vindicate a value of halo/disk normalization of 1/900, with probable root- 
mean-square range between 1/850 and 1/950, entirely consistent with more traditional 
determinations in the visible by Schmidt, and by Bahcall and colleagues. 

There are, of course, issues of uncertainty that strictly fall outside the context of 
SKY, such as a concern about the completeness of the complement of sources that SKY 
incorporates. Could there, for example, be a major population omitted from SKY that 
might dominate the faint NIR counts, such as red (more extreme than the SKY category 
of "M LATE V"), or brown dwarfs? While we cannot preclude this possibility, there are 
no indications of systematic under-predictions by SKY of observed counts at the faintest 
levels for which meaningful areas have so far been probed, namely K' <18 (of course, 
one must always be alert for the difficulties inherent in star-galaxy separation and for 
Malmquist bias in faint counts). Nor is there any evidence for the existence of any such 
cool, faint, halo population from studies with HST (Gould, Bahcall & Flynn 1996, 1997). 
Indeed, it might be more fruitful to invert the problem and to see which brown dwarf 
scenarios can already be precluded by using a model like SKY as an interpreter of exist- 
ing deep star counts. 

We thank the SSO staff for their technical support for these observations, and also thank 
K. Nakamura for her technical support of the observational software. YY acknowledges 

11 



financial support from the Yamada Science Foundation for the PICNIC camera and for 
the Japanese team's observing in Australia. TM was supported by the Grant in Aid 
for JSPS fellows by the Ministry of Education, Science, and Culture. This work has 
been supported in part by the Grant-in-Aid for Center-of-Excellence (COE) Research 
(07CE2002) of the Ministry of Education, Science, and Culture. The development of 
SKY has been supported largely by Dr. S. D. Price through the sponsorship of Phillips 
Laboratory under contract F19628-92-C-0090, with Vanguard Research, Inc. MC also 
thanks NASA- Ames for support under cooperative agreement NCC 2-142 with UC Berke- 
ley. 

References 

Arendt, R. G., Freundenreich, H. T., k Dwek, E. 1996, in "Unveiling the Cosmic Infrared 

Background", ed. E. Dwek, AIP Conf. Proc. 348, p.305 
Bahcall, J. N., Schmidt, M., k Soneira, R. M. 1983, ApJ, 265, 730 
Bahcall, J. N., & Soneira, R. M. 1984, ApJS, 55, 67 
Cohen, M. 1994, AJ, 107, 582 
Cohen, M. 1995, ApJ, 444, 874 

Cohen, M. 1997, in "The Impact of Large Scale Near-IR Sky Surveys", eds. F. Garzon et al. 

(Kluwer Academic), p. 67 
Cohen, M., Sasseen, T., k Bowyer, S. 1994, ApJ, 427, 848 

Cohen, M., Wainscoat, R. J., Walker, H. J., Volk, K., k Schwartz, D. E. 1989, AJ, 97, 1759 

Conti, P. k Vacca, W. 1990, AJ, 100, 431 

Gould, A., Bahcall, J.N. k Flynn, C. 1996, ApJ, 465, 759 

Gould, A., Bahcall, J.N. k Flynn, C. 1997, ApJ, 482, 913 

12 



Hammersley, P. L. 1996, priv. comm. 

Hammersley, P. L., Garzon, F., Mahoney, T., & Calbet, X. 1994, MNRAS, 269, 753 
Kobayashi, Y., Fang, G, Minezaki, T., Waseda, K. & Nakamura, K. 1994, Pore. SPIE, 2198, 
603 

McLeod, B.A., Bernstein, G.M., Rieke, M.J., Tollenstrup, E. V. & Fazio, G.G. 1995, ApJS, 
96, 117 

Meadows, V.S. 1994, Ph.D. dissertation, University of Sydney, Australia 

Minezaki, T., Kobayashi, Y., Yoshii, Y., & Peterson, B. A. 1997, accepted for ApJ (February 

10, 1998) 
Skrutskie, M. 1996, priv. comm. 
Schmidt, M. 1975, ApJ, 202, 22 

Wainscoat, R. J., Cohen, M., Volk, K., Walker, H. J., & Schwartz, D. E., 1992: ApJS, 83, 

111 [WCVWS] 
Yamagata, T. k Yoshii, Y. 1992, AJ, 103, 117 
Yoshii, Y. 1992, PASJ, 34, 365 

Yoshii, Y., Ishida, M, & Stobie, R. S. 1988, AJ, 93, 323 



13 



TABLE 1. The K' band star counts 
Kl Raw N a Completeness 6 n c Error c 

12.0-12.5 1 0.998 39.9 39.9 

12.5-13.0 1 0.998 39.9 39.9 

13.0-13.5 7 0.998 279. 106. 

13.5-14.0 6 0.998 240. 97.8 

14.0-14.5 3 0.998 120. 69.1 

14.5-15.0 5 0.998 200. 89.3 

15.0-15.5 11 0.998 439. 132. 

15.5-16.0 6 0.998 240. 97.8 

16.0-16.5 6 0.998 240. 97.8 

16.5-17.0 10 0.998 399. 126. 

17.0-17.5 11 0.995 440. 133. 

Notes to Table 1: 
a Raw counts of detected stars in the specified magnitude range. 
b The average of completeness for 15 < K' < 17 was presented at K' < 17. 
c Corrected star counts and the errors per magnitude per degree 2 . 
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Figure captions 



Figure 1: Our best match using SKY to the observed PICNIC star counts for 
halo/disk of 0.56 x the 1/500 value used by WCVWS, and z of 16.5 pc. Error bars are 
the actual Poisson errors associated with the observed star counts. The curves represent: 
the total (differential) counts (solid), and the separate disk (dotted, lighter curve below 
the total), and halo (long-dashed-dotted) counts. 

Figure 2: Quantitative estimates of the elements of fractional uncertainty associated 
with the extrapolation of our PICNIC corrected counts to faint K' levels. The "Total" 
curve is constructed from the root-sum-square of the 3 separate curves that arise from 
uncertainties in our ratio of halo/disk, in the space density of late-M dwarf stars, and in 
the location of the Sun with respect to the Galactic plane. 
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